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HIGHLIGHTS 


•  The  discharge  behavior  of  Li-S  cells  in  a  wide  range  of  current  rate  is  studied. 

•  Sulfur  must  highly  dissolve  in  electrolyte  to  maintain  high  capacity  performance. 

•  A  threshold  of  cathode’s  conductivity  exists  bellow  which  the  cell  would  not  operate. 

•  Above  the  conductivity  threshold  no  capacity  loss  is  observed  due  to  conductivity. 
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A  sensitivity  analysis  of  a  mathematical  model  of  a  lithium-sulfur  (Li— S)  battery  was  performed,  focusing 
on  the  effects  of  the  discharge  current  and  the  electronic  conductivity  of  the  cathode  by  investigating  the 
response  of  the  model  to  variations  of  these  parameters  over  a  wide  range  of  values.  The  behavior  of  the 
model  battery  system  with  respect  to  the  variations  of  each  parameter  is  explained  in  detail.  The 
nonlinear  response  of  the  model  to  the  parameter  variations  is  presented.  In  particular,  the  discharge 
current  causes  different  patterns  of  the  voltage  plateaus.  The  importance  of  the  electronic  conductivity  of 
the  cathode  is  also  discussed  in  detail.  A  minimum  conductivity  requirement  to  activate  the  cell  was 
observed,  beyond  which  no  significant  improvement  of  the  capacity  performance  occurs.  This  report 
evaluates  the  electrochemical  aspects  of  different  experimental  results  particularly  regarding  discharge 
voltage  plateaus,  and  moreover,  demonstrates  the  limits  of  the  model. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  lithium-sulfur  (Li-S)  battery  is  a  promising  energy  storage 
system  because  sulfur,  the  active  material,  has  a  high  theoretical 
specific  energy  (2600  Wh  kg-1);  furthermore,  sulfur  is  inexpensive, 
abundant,  and  nontoxic.  Over  the  past  decade,  extensive  efforts 
have  focused  on  developing  a  rechargeable  Li— S  battery  [1].  How¬ 
ever,  no  Li-S  battery  has  been  commercialized  to  date  due  to 
several  unsolved  problems.  Because  of  the  insulating  nature  of  the 
Li-S  battery  discharge  products,  namely,  sulfur  and  lithium  sulfides 
(Li2S2,  Li2S),  a  relatively  high  amount  of  conductive  material  must 
be  added  to  the  Li-S  battery  cathode  to  guarantee  the  utilization  of 
the  active  material  [2].  In  addition,  the  dissolved  polysulfides  in  the 
electrolyte  diffuse  through  the  separator  to  the  lithium  anode  and 
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react  directly  with  the  lithium.  This  reaction  creates  the  so-called 
internal  shuttle  phenomenon,  which  causes  passivation  of  active 
material  and  the  self-discharge  of  the  battery  [1-3]. 

Although  a  wide  range  of  experimental  strategies  have 
attempted  to  tackle  these  problems,  little  theoretical  attention  has 
been  paid  to  understanding  the  mechanism  of  the  Li— S  battery. 
Moreover,  the  behavior  of  the  sulfur  as  the  active  material  is  very 
complicated  and  extremely  sensitive  to  the  different  physical  and 
chemical  parameters  involved  in  the  preparation  of  the  cathode 
composite  or  even  during  the  assembly  of  a  cell  [2,4-6].  For 
example,  very  complex  reduction  reactions  of  elemental  sulfur,  Sg, 
occur  during  the  discharge  processes,  and  different  polysulfide 
chains  form  in  as-yet-undetermined  ways  [4-6].  Although  the 
exact  reduction  reaction  mechanism  is  not  yet  definitively  deter¬ 
mined,  some  models  have  been  introduced  to  describe  Li-S  bat¬ 
teries.  Y.  V.  Mikhaylik  et  al.  [7]  developed  a  mathematical  model  for 
the  shuttle  effect  based  on  a  two-stage  reduction  of  sulfur.  Specif¬ 
ically,  Sg  reduces  to  S4-,  and  S4-  reduces  to  S2-  and  S2-.  The  direct 
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Fig.  1.  Schematic  of  a  Li-S  cell  and  a  summary  of  the  governing  equations  and  boundary  conditions. 


reduction  of  high-polysulfides  on  the  surface  of  the  lithium  anode 
was  also  included  in  the  model.  The  model  produced  a  good  un¬ 
derstanding  of  the  shuttle  effect,  thus  providing  evidence  that  the 
self-discharge,  overcharge,  and  efficiency  are  strongly  related  to  the 
shuttle  effect. 

Another  mathematical  model,  introduced  by  K.  Kumaresan  et  al. 
[8],  considered  five  different  stages  for  the  reduction  of  sulfur  to 
mathematically  describe  the  discharge  behavior  of  Li-S  batteries. 
The  model  considers  the  dissolution  and  precipitation  of  sulfur  and 
polysulfides  into  the  electrolyte.  Assuming  that  the  electrochemical 
reactions  in  the  system  follow  the  Butler— Volmer  equation,  the 
equations  governing  the  material  balance  of  each  individual  species 
are  solved  in  the  model.  The  authors  assumed  many  different  pa¬ 
rameters  that  must  be  determined  by  proper  experiments,  but  in 
the  absence  of  such  experimental  results,  they  provided  educated 
guesses  for  most  of  the  parameters.  Nevertheless,  their  results 
corresponded  well  with  a  class  of  experimental  results. 

In  fact,  the  behavior  of  sulfur  batteries  strongly  depends  on  the 
morphology  [9]  and  chemical  properties  of  the  cathode  composite 
and  electrolyte  chosen  [2].  Therefore,  changing  the  additives  in  the 
cathode  composite  [2,9,10]  or  electrolyte  [2,11]  leads  to  different 
discharge  plateau  shapes.  Two  flat  discharge  plateaus  are  observed 
in  most  cathode  composites  of  sulfur  and  carbon  materials  [2  .  In 
contrast,  in  the  cathode  composites  made  of  sulfur  and  polymer 
materials,  only  one  decreasing  plateau  typically  appears  in  the 
discharge  duo  in  the  interaction  between  sulfur  and  polymer 
[10,12].  Even  within  the  two  categories,  significant  changes  occur  in 
the  details  of  the  discharge  plateau  in  different  experiments. 
Despite  the  difficulty  of  doing  so,  investigating  a  specific  system 
using  the  model  of  Kumaresan  et  al.  requires  that  each  parameter 
be  determined  through  an  analysis  of  the  results  of  experiments  on 
that  system.  Conversely,  because  of  the  system’s  complexity, 
different  sets  of  parameters  may  result  in  the  same  model  output.  A 
sensitivity  analysis  of  the  model  parameters  will  help  elucidate  the 
model  and  identify  the  actual  Li-S  cell  behavior.  The  results  of  a 
sensitivity  analysis  would  identify  a  range  of  values  for  each 
physical  parameter  and,  more  importantly,  a  set  of  possible  physical 
or  chemical  properties  responsible  for  a  specific  behavior  of  the  cell. 

The  aim  of  this  paper  is  not  only  to  investigate  the  possible 
ranges  of  the  physical  parameters  by  performing  a  sensitivity 


analysis  but  also  to  explore  the  mathematical  validation  of  the 
model  and  thus  clarify  its  abilities  and  limits.  The  results  of  the 
sensitivity  analysis  can  be  used  to  improve  the  mathematical  model 
and  can  ultimately  contribute  to  the  development  of  commercial 
Li— S  batteries. 

This  article  is  the  first  part  of  a  series  on  the  sensitivity  analysis 
of  the  model  presented  by  Kumaresan  et  al.  In  this  article,  the  focus 
is  on  the  response  of  the  model  to  variations  in  the  discharge  cur¬ 
rent.  Furthermore,  the  electronic  conductivity  of  the  cathode,  one 
of  the  major  issues  in  Li-S  batteries,  is  studied  in  detail.  The  dy¬ 
namics  of  the  active  material  in  the  solid  phase  is  investigated  to 
estimate  the  cyclability  of  the  Li-S  cell. 

2.  Model  development  and  governing  equations 

The  model  is  based  on  the  proposed  model  by  K.  Kumaresan,  Y. 
Mikhaylik,  and  R.  E.  White  [8].  Fig.  1  presents  a  schematic  of  the  Li- 
S  cell  and  summarizes  the  governing  equations  and  boundary 
conditions  of  the  model.  Appendix  B  describes  the  model,  and 
Appendix  A  defines  its  parameters.  The  values  proposed  for  the 
parameters  are  provided  in  Tables  1  -4  in  Appendix  C. 

The  reactions  in  the  cell  include  Li  metal  oxidization  at  the 
anode  surface  during  discharge  [8]: 

Li^Li++e-  (1) 

It  is  assumed  that  during  discharge,  the  elemental  sulfur,  which 

is  initially  in  the  solid  phase,  dissolves  in  the  electrolyte  and  then 

goes  through  the  following  electrochemical  reactions  8]: 

ls8(1)+e—  ls|-  (2) 

|s|-+e^2S|-  (3) 

Sl+e-^sj-  (4) 

\s2A-+e-^S22-  (5) 
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^sj-+e-^±2S2-  (6) 

The  following  precipitation/dissolution  reactions  are  also 


present: 

58(5)^58(1) 
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The  details  of  the  model  and  its  governing  equations  can  be 
found  either  in  Appendix  B  or  Ref.  [8].  The  current  densities  due  to 
the  electrochemical  reactions  are  given  by  the  Butler— Volmer 
equation. 

3.  Results  and  discussion 

A  range  of  variations  is  assumed  for  each  parameter.  This  range 
is  not  bounded  by  a  range  of  physical  values.  The  goal  is  to  deter¬ 
mine  the  behavior  of  the  model  system  with  respect  to  different 
situations  and  also  to  find  a  range  for  the  parameters  in  which  the 
Li-S  cells  are  feasible.  However,  the  functionality  of  the  model  with 
respect  to  these  parameters  was  found  not  to  be  linear.  Instead,  we 
must  choose  a  parameter  each  time  and  investigate  the  changes  in 
the  model  from  variations  in  this  parameter  while  keeping  the 
other  parameters  constant. 

3.1.  Applied  discharge  current ,  Iapp 

One  of  the  important  goals  in  battery  research  is  to  widen  the 
discharge  current  operation  range  of  batteries;  for  example,  a 
particularly  high  current  rate  is  required  for  electric  vehicles. 
Moreover,  the  discharge  voltage  plateau  of  a  battery  and  its  cycla- 
bility  are  strongly  dependent  on  the  current  rate  during  operation. 
Therefore,  the  response  of  the  model  to  a  range  of  applied  currents 
must  be  investigated. 

Fig.  2  presents  the  discharge  voltage  plateaus.  In  the  first 
three  discharge  currents,  with  low  applied  discharge  currents, 
Japp  =  0.02C,  0.05C,  and  0.1  C,  two  voltage  plateaus  are  very  clear, 
and  the  voltage  curves  are  quite  similar,  except  from  a  decreasing 
shift  in  potential.  This  similarity  indicates  that  the  kinetics  of  the 


Fig.  2.  The  discharge  plateau  of  a  sulfur-based  cell  at  different  discharge  current  rates. 
First  plateau  vanishes  due  to  low  dissolution  rate  of  elemental  sulfur. 


chemical  reactions  is  almost  the  same  at  these  rates  of  discharge. 
Because  the  decrease  in  the  potential  due  to  the  resistivity  of  the 
cathode  matrix  and  electrolyte  is  very  small  (<1(T4  V)  compared  to 
the  total  decrease  in  potential,  one  may  conclude  that  the  kinetics 
of  the  chemical  reactions  causes  the  decrease.  In  other  words,  the 
decrease  is  due  to  polarization  at  the  surface  of  the  anode  and 
cathode.  The  linear  approximation  of  the  Butler— Volmer  equation 
(Equation  (B5))  describes  a  linear  relationship  between  the 
decrease  in  potential  and  change  in  current  (Ai°c  AV)  and  thus  well 
approximates  the  observed  plateaus  in  the  first  three  discharge 
currents. 

The  first  plateau  starts  to  deform  at  the  discharge  rate  of  0.05C, 
as  the  speed  of  dissolution  of  elemental  sulfur  into  the  electrolyte  is 
too  slow  to  maintain  a  constant  concentration  of  dissolved  sulfur. 
The  slow  dissolving  rate  of  sulfur  is  the  main  reason  for  the  capacity 
loss,  which  starts  at  2C.  Fig.  3  provides  more  details  of  the  reactions 
kinetics  of  the  Li-S  system.  The  normalized  current,  which  is  due  to 
electrochemical  reaction  j  in  the  cathode  and  is  defined  by  the 
following  expression,  is  used  to  illustrate  the  system’s  kinetics: 

x  =  L 

If  =  A-  [  aijdx  (12) 

fapp  J 
x  =  Ls 

where  Y%=  2$  =  1.  In  other  words,  1 rf*  represents  the  overall 
contribution  of  the  electrochemical  reaction  j  in  producing  the 
discharge  current  at  a  given  time.  In  the  ideal  case  of  utilizing  the 
entire  capacity  of  sulfur,  the  discharge  normalized  currents  satisfy 
the  following  relations: 

Tdc 

2-  [  Ifdt  =  Cj,  (13a) 

JDC  i  J 
t  =  0 

<c* . «  -  {wi«)  (,3b> 

where  Tdc  is  the  total  time  of  an  ideal  complete  discharge  and  Cj  is 
the  portion  of  reaction  j  in  the  total  capacity.  The  average  concen¬ 
tration  of  the  sulfide  species  at  the  cathode  is  presented  in  Fig.  4. 
Because  the  graphs  of  normalized  current  for  the  cases  Japp  =  0.1  C 
and  0.05C  are  similar  to  the  case  of  0.02C,  they  are  not  shown  in 
Fig.  3,  whereas  the  concentration  of  the  sulfide  species  in  these 
cases  are  shown  (dotted  lines)  to  illustrate  the  difference  in  the 
second  part  of  the  discharge.  As  the  discharge  current  rate  increases 
to  0.1  C,  the  concentration  of  the  high  sulfides  (S8(i>  S|+,  and  Sg+) 
decreases  more  in  the  second  part  than  in  the  first  part  because 
more  material  is  used  in  the  reaction.  Thus,  a  higher  concentration 
of  the  low  sulfides  (S|+,  S2+,  and  S2+)  forms  before  the  rate  of  the 
low-sulfide  reactions  becomes  sufficiently  high  and  the  precipita¬ 
tion  of  Li2S  begins.  This  sequence  of  events  causes  a  small  shift  in 
the  sharp  trough  between  the  two  stages.  The  shift  in  precipitation 
toward  the  right  occurs  even  for  a  discharge  current  rate  of  0.5C,  for 
which  the  trough  starts  to  shift  because  the  low-sulfide  reaction 
begins  relatively  sooner.  The  relatively  fast  resolution  of  sulfur  to 
electrolyte  maintains  a  constant  sulfur  concentration  during  the 
first  stage,  and  therefore,  a  constant  reaction  rate  was  observed.  The 
reduction  of  occurs  simultaneously  with  a  constant  rate 
because  the  reference  voltage  is  similar.  When  the  solid  sulfur  in 
the  system  is  fully  consumed,  the  concentration  of  sulfur  decreases 
suddenly  (in  the  related  time  scale),  as  does  the  related  current.  At 
this  point,  the  second  reaction  reaches  a  sharp  peak  in  its  time  scale 
and  quickly  reduces  almost  all  of  the  S .  The  reduction  of  Sg_ 
becomes  the  dominant  reaction  in  the  system,  leading  to  the 
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Fig.  3.  (a),  (b),  (c),  (d),  and  (e)  Normalized  current  (see  Equation  (12)),  due  to  different  electrochemical  reactions,  produced  by  the  reduction  of  S8(i),  S|_,  S§-.  S4-,  and  S^-, 
respectively,  (f)  Average  volume  fraction  of  solid  elemental  sulfur,  es  and  Li2S(S),  £UaS  ,  in  the  cathode  during  discharge.  Unutilized  elemental  sulfur  is  the  main  reason  of  capacity 
lost  in  high  C  rates. 


production  of  the  lower  sulfides.  As  a  result,  the  concentration  of 
the  low  sulfides  increases;  therefore,  the  reaction  requires  a  higher 
polarization  at  the  cathode’s  surface.  The  potential  of  the  cell 
consequently  decreases  to  a  local  minimum  in  a  sharp  trough.  At 
the  point  at  which  Li2S  starts  to  precipitate,  the  concentration  of 
S2-  decreases,  causing  a  sharp  increase  in  rate  of  the  last  two  re¬ 
actions.  Interestingly,  because  the  S4-  concentration  is  high  for  the 
related  potential,  a  reverse  reaction  occurs  for  the  high  sulfides. 
Upon  reaching  a  quasi-steady  state,  the  rate  of  the  last  two  re¬ 
actions  becomes  constant  during  the  remaining  discharge  time. 
However,  the  third  reaction  continues  slowly.  Once  all  of  the  sul¬ 
fides  have  been  reduced,  a  sudden  drop  in  the  potential  occurs. 

For  a  discharge  current  of  0.5C  and  higher,  sulfur  dissolution  is 
too  slow  to  maintain  a  constant  concentration  of  elemental  sulfur 
in  the  electrolyte  solution.  Thus,  the  other  reactions  begin  corre¬ 
spondingly  sooner  in  their  time  scales,  creating  smoother  peaks 
and  drops.  However,  up  to  a  1C  discharge  rate,  the  dissolution  is 
sufficiently  fast  to  keep  the  first  (and  second)  reaction  dominant 
during  the  first  stage  of  discharge,  but  those  reactions  continue  in 
the  second  part  of  the  discharge  because  some  solid  sulfur  is  still 
present  that  can  be  dissolved  into  the  electrolyte  solution.  The 
reverse  reactions  then  cease  due  to  the  earlier  start  of  the  low- 
sulfide  reactions  in  their  relative  time  scales.  Although  the  sharp 
troughs  in  potential  and  in  the  concentrations  of  high  sulfides  and 
concentration  peaks  of  the  low  sulfides  still  exist  when  S2~  pre¬ 
cipitation  starts  (for  the  reasons  explained  earlier),  the  sharp 
trough  in  the  high-sulfide  concentrations  is  no  longer  due  to  a 


reverse  reaction.  Instead,  this  trough  now  occurs  as  the  species 
consumption  rate  drops  due  to  a  reduction  reaction,  concurrent 
with  the  production  due  to  the  dissolving  of  sulfur.  The  dependence 
of  the  dissolution/precipitation  of  each  species  on  the  nucleus  size 
of  each  precipitate  is  clearly  observed  in  the  volume  fraction  of 
elemental  sulfur,  which  does  not  exhibit  a  linear  decrease. 

Starting  from  2C,  elemental  sulfur  does  not  dissolve  sufficiently 
fast  to  maintain  even  the  first  two  reactions  as  the  system’s 
dominant  reactions.  The  low-sulfide  reductions  start  earlier,  and 
the  low-sulfide  concentrations  reach  their  maximum  levels  very 
quickly.  Again,  the  start  of  the  S2~  precipitation  causes  the  sharp 
trough  in  the  levels  of  the  high-sulfide  concentrations  and  the 
peaks  in  low-sulfide  concentrations.  Whereas  the  high-sulfide  re¬ 
action  rates  change  smoothly  over  time,  the  reduction  of  S4-  ex¬ 
hibits  a  sharp  trough  at  this  starting  point  of  precipitation  and  the 
last  reaction  reaches  a  peak.  The  dissolution  of  elemental  sulfur 
does  not  occur  sufficiently  quickly:  the  solid  sulfur  particles 
dissolve  more  slowly  as  they  become  smaller.  Thus,  the  sulfide 
species  cannot  be  produced  in  the  quantities  required  to  be 
involved  in  the  reactions,  causing  the  potential  to  drop  and  the  end 
of  discharge,  and  some  unused  solid  sulfur  remains  in  the  cathode. 

At  low  discharge  rates,  the  species  concentration  is  nearly  uni¬ 
form  in  the  cell  because  diffusion  proceeds  faster  than  the  pro¬ 
duction  and  consumption.  However,  at  high  C  rates,  gradients  in  the 
concentration  of  the  species  are  generated  across  the  cell  to  pro¬ 
duce  the  essential  mass  transfer.  For  all  the  species,  except  S2-  and 
Sg(i),  the  maximum  concentration  is  at  the  interface  of  the  anode 
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- 0.02  C  - 0.5  C  - 1C  - -2  C  - 3  C  - 5C  - 7C  0.05  C  0.1  C 


Fig.  4.  Average  concentration  of  the  relevant  species  in  the  cathode  during  discharge.  Increasing  the  discharge  rate,  the  concentrations  of  high  polysulfides  decrease  and  the 
concentrations  of  low  polysulfides  increase. 


and  separator.  The  concentration  gradients  have  higher  slopes  at 
the  separator  than  at  the  cathode.  The  S2_  concentration  exhibits 
the  same  behavior  before  precipitation  starts  but  drops  at  both  the 
separator  and  cathode  once  precipitation  starts.  However, 
comparing  the  concentrations  in  both  areas,  the  S2-  concentration 
remains  higher  at  the  cathode  because  of  its  simultaneous  pro¬ 
duction.  Consequently,  a  concentration  gradient  forms  at  the 
cathode-separator  interface.  In  contrast,  a  high  concentration  of  S2- 
and  Li+  at  the  anode  and  separator  interface  causes  a  relatively  high 
precipitation  of  Li2S(S)  at  the  anode  surface  soon  after  discharge 
begins.  Fig.  5  shows  the  volume  fraction  of  Li2S(S)  across  the  cell  at 
different  discharge  times,  at  Japp  =  7C.  At  t  =  10  s,  the  volume 
fraction  is  almost  zero  everywhere.  At  t  =  20  s,  the  volume  fraction 
begins  to  increase  and  has  a  large  maximum  at  the  anode  surface. 
Over  time,  the  precipitation  continues  everywhere  in  the  cathode 
because  the  production  of  S2-  occurs  at  that  location.  Moreover,  the 
diffusion  of  S2-  into  the  separator  causes  Li2S(S)  to  form  deep  in  the 
separator  (for  the  same  reason,  the  volume  fraction  in  the  cathode 
at  a  short  distance  from  the  interface  is  less  than  that  deep  inside 
the  cathode).  However,  the  higher  concentration  of  Li+  at  the 
interface  causes  a  Li2S(S)  volume  fraction  peak  to  form  at  a  short 
distance  from  the  interface  inside  the  cathode.  This  information 
elucidates  the  stresses  in  the  cell  and  the  deformation  of  the 
cathode  during  cycling. 

3.2.  Conductivity  of  the  cathode,  a 

Neither  sulfur  nor  polysulfides  are  conductive;  therefore,  the 
cathode  must  consist  of  a  conductive  porous  matrix  that  holds 


sulfur  in  its  pores  and  facilities  the  accessibility  of  electrons  to  the 
sulfur.  The  method  of  manufacturing  the  conductive  matrix  is  one 
of  the  major  difficulties  in  sulfur-based  cell  fabrication:  on  the  one 
hand,  a  higher  conductivity  requires  more  conductive  additives  and 
less  sulfur  in  the  cathode,  causing  the  cell  to  lose  its  specific  energy 
capacity;  on  the  other  hand,  lower  conductivity  causes  reduced 
utilization  of  the  sulfur  as  the  active  material.  This  section  in¬ 
vestigates  the  effect  of  the  conductivity  of  the  porous  matrix  on  the 
cell  performance.  No  difference  in  the  matrix  porosity  is  assumed; 
only  the  conductivity  was  varied.  This  investigation  approach 
matches  the  situations  in  which  the  conductive  additive  is  changed 
but  not  the  amount  used,  assuming  the  porosity  remains  the  same. 


Gpp  —  7  C 


Fig.  5.  Volume  fraction  of  Li2S(S),  £Li2s(s)  across  the  cell  at  different  times  of  discharge 
with  a  current  rate  of  7C.  The  interface  of  the  separator  and  the  cathode  is  located  at 
x  =  9  pm. 
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Fig.  6  illustrates  the  effect  of  the  conductivity  of  the  cathode 
porous  media  matrix  on  the  battery  performance.  The  upper  and 
lower  sets  of  lines  represent  the  voltage  plateaus  for  a  1C  and  5C 
discharge  rate,  respectively.  The  conductivity  of  the  cathode 
matrix  is  assumed  to  vary  from  1  S  m-1  to  10-4  S  m-1.  For  the 
discharge  rate  of  5C,  the  simulation  result  for  a  =  10-4  S  m-1  did 
not  converge  to  a  solution  at  times  longer  than  approximately 
2  s,  i.e.,  it  is  not  possible  to  discharge  at  a  high  C  rate  when  the 
conductivity  is  not  sufficiently  high.  The  graphs  indicate  that 
once  a  certain  conductivity  threshold  has  been  exceeded,  no 
significant  change  occurs  to  the  voltage  plateau  if  the  conduc¬ 
tivity  is  increased  further,  indicating  the  dominant  role  of  the 
reaction  kinetics  in  the  cell  behavior.  As  the  conductivity  de¬ 
creases,  an  unexpected  drop  in  potential  at  the  beginning  of  the 
cell  discharge  was  observed.  This  drop  is  more  obvious  at  a  rate 
of  1C.  At  low  conductivity,  the  potential  gradient  in  the  cathode 
matrix  is  steep  (unlike  at  high  conductivity),  with  a  minimum  at 
the  interface  of  the  cathode  and  current  collector.  Therefore,  a 
steep  gradient  in  the  concentration  of  the  species  is  required  for 
a  well-distributed  electrochemical  reaction  across  the  cathode. 
As  the  cell  discharge  begins  because  the  species  concentration  is 
uniform,  a  sharp  drop  in  the  potential  of  the  solid  matrix  forms  at 
the  interface  of  the  cathode  and  current  collector  to  draw  all  of 
the  current  from  this  region.  Over  time,  the  production  of  more 
ions  in  the  electrolyte  causes  a  gradient  in  the  concentration  of 
species  to  form,  leading  to  more-distributed  electrochemical  re¬ 
actions  across  the  cathode.  However,  the  peak  of  the  reactions 
still  exists  at  that  interface,  but  it  is  not  as  sharp  as  before.  At  the 
end  of  the  discharge,  the  smaller  number  of  ions  again  causes  to 
a  sharp  potential  drop  to  form  at  the  interface.  The  high  rate  of 
low-sulfide  production  at  the  interface  leads  to  a  high-rate  of 
precipitation  as  well  such  that  precipitated  species  fill  the  porous 
regions.  Fig.  7  presents  the  volume  fraction  of  Li2S(S)  at  the 
separator  and  cathode  during  discharge.  The  porous  medium,  at 
the  interface  of  the  cathode  and  current  collector,  is  clearly 
completely  filled  by  the  precipitant.  One  may  conclude  that  this 
filling  phenomenon  causes  the  active  material  to  be  lost  in  the 
cycling  because  the  solid  sulfur  at  that  area  does  not  have  easy 
access  to  the  electrolytes. 

4.  Conclusions 

A  sensitivity  analysis  of  a  mathematical  model  of  a  Li-S  cell  with 
respect  to  the  discharge  current  rate  and  the  conductivity  of  the 
cathode  matrix  was  performed,  the  results  of  which  provided  de¬ 
tails  on  the  different  characteristics  of  the  cell.  Depending  on  the 
discharge  current  (and  the  availability  of  dissolved  sulfur),  the 
relevant  electrochemical  reactions  can  occur  either  simultaneously 


Fig.  6.  The  discharge  plateaus  of  the  sulfur-based  cell  with  different  conductivities  of 
the  cathode  matrix  for  discharge  current  rates  of  5C  and  1C. 


Ippp  =  5  C;  c=  W3  S/m 


Fig.  7.  Volume  fraction  of  Li2S(S),  cU2S(s)  across  the  cell  at  different  times  of  discharge  for 
a  discharge  current  rate  of  5C.  The  separator-cathode  interface  is  located  at  x  =  9  pm. 


or  after  the  previous  reactions  have  completed.  The  coincidence  or 
non-coincidence  of  the  electrochemical  reactions  determines  the 
shape  of  the  discharge  plateau.  Sharp  changes  in  the  voltage 
plateau  are  observed  when  the  electrochemical  reactions  occur 
non-coincidentally,  whereas  a  smooth  plateau  is  expected  in  the 
coincidental  case.  In  particular,  if  the  first  plateau  has  vanished,  one 
may  verify  the  existence  of  the  first  plateau  by  decreasing  the 
discharge  current.  If  the  plateau  appears  at  a  low  discharge  rate, 
slow  dissolution  of  elemental  sulfur  into  the  electrolyte  causes  a 
significant  capacity  loss.  The  details  of  the  effects  of  precipitation 
rate  constant  on  the  cell  behavior  will  be  studied  in  part  III  of  this 
series  of  papers.  The  simulation  results  also  demonstrate  a  partial 
movement  of  the  active  material  into  the  separator  during 
discharge. 

To  functionalize  the  cell,  a  minimum  in  the  cathode  conductivity 
is  required  (depending  on  the  discharge  current)  below  which  the 
cell  would  not  operate  and  above  which  no  capacity  loss  is 
observed  due  to  the  conductivity.  Low  conductivity  only  causes  a 
steeper  discharge  voltage,  particularly  at  the  beginning  and  end  of 
discharge. 

The  model  predicts  a  very  flat  second  plateau,  which  differs 
from  the  results  of  many  experiments.  This  difference  is  a  direct 
consequence  of  assuming  that  all  the  surfaces  in  the  cathode  are 
active  surfaces  for  electrochemical  reactions,  ignoring  the  non- 
conductive  nature  of  the  solid  sulfur  and  polysulfides.  The  model 
needs  to  be  improved  by  including  active  surface  loss  due  to  the 
precipitation  of  polysulfides. 
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Appendices 

A.  Parameters  and  symbols  (primarily  taken  from  Ref  [5]) 


a  specific  surface  area  of  the  cathode 

a0  initial  value  of  a 

b  Bruggeman  coefficient 

Q  concentration  of  species  i 

(i  =  Li+ ,  S8(</) ,  S|— ,  S§— ,  ,  S|— ,  S2— ,  and  A~  (anion  of  the 

lithium  salt  used  in  the  electrolyte)),  mol  itt3 
Q.ref  reference  concentration  of  species  i,  mol  m-3 
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Cj  portion  of  reaction  j  in  the  total  capacity  for  an  ideal 

complete  discharge 

D/  o  diffusion  coefficient  of  species  i  in  the  bulk  medium, 

2  — 1 

m  s 

D  diffusion  coefficient  of  species  i  in  the  porous  medium, 

2  — 1 

m  s 

F  Faraday  constant,  C  equi-1 

ij  current  density  due  to  reaction  j,  A  m-2 

ij>ref  exchange  current  density  of  the  electrochemical  reaction  j 

at  the  reference  concentrations,  A  m-2 
i\  superficial  current  density  in  the  liquid  phase,  A  nrr2 

is  superficial  current  density  in  the  solid  phase,  A  m-2 

Japp  applied  current  density,  A  m-2 
1^  normalized  current  due  to  electrochemical  reaction  j 

I< sPfk  solubility  product  of  precipitate  k 

ki<  rate  constant  of  precipitate  k 

L  thickness  of  the  cell,  m 

Ls  thickness  of  the  separator,  m 

Nj  superficial  flux  of  species  i,  mol  m2  s_1 

tij  number  of  electrons  transferred  in  electrochemical 

reaction  j 

Pij  anodic  reaction  order  of  species  i  in  electrochemical 
reaction  j 

qij  anodic  reaction  order  of  species  i  in  electrochemical 
reaction  j 

R  gas  constants  mol-1  K-1 

Ri  production  rate  of  species  i  due  to  precipitation  reactions, 

mol  m3  s-1 

R'k  rate  of  precipitation  of  solid  species  k,  mol  m3  s-1 
r*  production  rate  of  species  i  due  to  electrochemical 

reactions,  mol  m3  s_1 

Sij  stoichiometric  coefficient  of  species  i  in  electrochemical 
reaction  j 

T  temperature,  K 

Tdc  total  time  of  an  ideal  complete  discharge,  s 
t  time,  s 

Uj  standard  open  circuit  potential  (OCP)  of  electrochemical 

reaction  j 

Ujt ref  OCP  of  electrochemical  reaction  j  at  reference 
concentrations,  V 

Vk  Molar  volume  of  the  precipitate  k,  m3  mol-1 
z/  charge  number  of  species  i 

aaj  anodic  transfer  coefficient  of  reaction  j 

aCj  cathodic  transfer  coefficient  of  reaction  j 
e  porosity  of  the  separator  and  cathode 

£k  volume  fraction  of  precipitate  k  in  the  separator  and 

cathode 

cp\  potential  in  the  liquid  phase,  V 

(ps  potential  in  the  solid  phase,  V 

jitk  number  of  ionic  species  i  produced  by  dissociation  of 
precipitate  k 

r\j  overpotential  for  electrochemical  reaction  j 

a  effective  conductivity  of  the  solid  phase  of  the  cathode, 

S  rrT1 

?  morphology  parameter 

B.  Governing  equations 


where  the  flux  is  given  by 

y  =  -D,VC,  -  z,^FC,VVl  (B2) 

The  diffusion  coefficient  for  species  i,  Du  is  corrected  based  on 
Bruggeman’s  expression  for  porosity  and  tortuosity:  Df  =  Di0£b. 
The  rate  of  production/consumption  of  species  i  due  to  the  elec¬ 
trochemical  reactions  can  be  written  in  the  form 


ri  =  aJ2S-4 

^  n.F 


(B3) 


where  the  stoichiometric  coefficients  are  given  in  Table  1,  Appendix 
C.  The  specific  surface  area  of  the  cathode  varies  due  to  the  pre¬ 
cipitation/dissolution  of  the  various  lithium  sulfide  species  and  is 
assumed  to  be  governed  by  the  empirical  expression: 


a  =  a0 


^initial 


(B4) 


where  the  empirical  parameter,  ?,  is  assigned  a  value  of  1.5.  The 
Butler-Volmer  equation  yields  the  current  density  due  to  the 
electrochemical  reactions: 


(B5) 


where  the  overpotential  for  the  reaction  j  is 


Vj  =  <Ps  -  <P\  -  Vj, ref  (B6) 

The  terms  =  Sy  refer  to  anodic  species  and  =  -Sy  refer  to 
cathodic  species.  The  open-circuit  potential  for  reaction  j  is  given  by 


RT 


"M,  Uf  n.F^S(. 


In 


Q,ref 

1000 


The  liquid  phase  current  density  is  given  by 


(B7) 


<1  =  FjfZiNi  (B8) 

i 

The  solid  phase  current  density  follows  Ohm’s  law 

is  =  -avcps  (B9) 

Because  the  charge  can  enter  or  leave  the  liquid  phase  only  by 
electrochemical  reactions,  the  following  equations  apply  at  the 
liquid/solid  interface: 


V'l  =  aJ2‘j 

j 


(B10) 


The  model  includes  the  electrochemical  reactions  given  by 
Equations  (1)— (6)  and  the  precipitation/dissolution  reactions  given 
by  Equations  (7)— (11 ).  In  a  porous  medium,  the  governing  equation 
for  the  material  balance  of  an  individual  species  is 

^p=*  -V  Ni  +  rt-Ri  (Bl) 


and 

V  •  is  +  V  •  i\  =  0  (Bll) 

The  rate  of  consumption  or  production  of  species  i  due  to  pre¬ 
cipitation/dissolution  is  related  to  the  rate  of  precipitation/disso¬ 
lution  reaction  k  by 
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Ri  =  Erf  (B12) 

k 

where  the  rate  of  precipitation  of  solid  species  k  is  assumed  to  be 
governed  by  the  following  kinetic  equation: 

(Bi3) 

The  volume  fraction  of  the  precipitate,  k,  is  a  function  of  time 


and  is  given  by 

d£k  \r  n/ 

IF  VkRk 

(B14) 

Therefore,  the  porosity  variation  with  time  is 

5-?W 

(B15) 

The  boundary  conditions  of  the  model  are  as 
interface  of  the  cathode  and  the  current  collector,  x 

follows.  At  the 
=  L,  the  BCs  are 

Nt  =  0 

(B16a) 

Is  =  Japp 

(B16b) 

h  =  o 

(B16c) 

At  the  cathode-separator  interface,  x  =  Ls, 

Nf, seperator  =  ^i, cathode 

(B17a) 

h,  seperator  =  h,  cathode  =  Japp 

(B17b) 

is  =  0 

(B17c) 

At  the  surface  of  anode,  x  =  0, 

(ps  =  0 

(B18a) 

Nt  =  0 

(B18b) 

Ni  =  h/F 

(B18c) 

ii  =  FN ! 

(B18d) 

The  governing  equations  can  be  solved  numerically. 
C.  Parameter  tables  for  reference  values 


Table  1 

Stoichiometric  coefficients,  Sy  [8]. 


Sy  Reactions  given  by  equation  ( j ) 


Species  (i) 

1 

2 

3 

4 

5 

6 

Li+ 

-1 

0 

0 

0 

0 

0 

S8(l) 

0 

-1/2 

0 

0 

0 

0 

s|- 

0 

1/2 

-3/2 

0 

0 

0 

1" 

s24- 

0 

0 

2 

-1 

0 

0 

0 

0 

0 

3/2 

-1/2 

0 

c2— 

0 

0 

0 

0 

1 

-1/2 

s2- 

0 

0 

0 

0 

0 

1 

A~ 

0 

0 

0 

0 

0 

0 

Table  2 

Kinetic  and  thermodynamic  properties  [8]. 


Reaction  (j) 

'ojref  <Am~2) 

<j 

nj 

uf 

1 

0.394 

0.5 

0.5 

1 

0.0 

2 

1.9719 

0.5 

0.5 

1 

2.39 

3 

0.019719 

0.5 

0.5 

1 

2.37 

4 

0.019719 

0.5 

0.5 

1 

2.24 

5 

1.972  x  10"4 

0.5 

0.5 

1 

2.04 

6 

1.972  x  10"7 

0.5 

0.5 

1 

2.01 

Table  3 

Transport  properties  and  reference  concentrations  [8]. 

Species  (i) 

Zi 

D*i0( m2  s-1) 

Cr  ref  (mo1  m-3) 

Li+ 

+1 

1  x  10"10 

1001.08 

S8(i) 

0 

10  x  10  10 

19.0 

sj- 

-2 

6  x  10"10 

0.1832 

Sl 

-2 

6  x  lO"10 

0.3351 

Sl 

-2 

1  x  lO"10 

0.02146 

T 

S2" 

-2 

1  x  lO"10 

5.999  x  10~7 

-2 

1  x  lO"10 

9.94  x  lO"10 

A~ 

-1 

4  x  lO"10 

1000 

Table  4 

Separator  and  cathode  parameters  [8]. 

Parameter 

Separator 

Cathode 

Thickness  (m) 

9  x  10“6 

41  x  10"6 

^initial 

0.37 

0.778 

£S8(s)  Initial 

1  x  10"12 

0.160 

£Li2S8(S), initial 

1  x  10“6 

1  x  10"6 

£Li2S4(S),  initial 

1  x  10“6 

1  x  10"6 

£Li2S2(s)  Initial 

1  x  10"6 

1  x  10"6 

£Li2S(S)  Initial 

1  x  lO"7 

1  x  10"7 

a0 

- 

132,762 

Table  5 

Parameters  for  the  precipitation  reactions,  either  assumed  or  from  Ref.  [8]. 


Precipitate 

(k) 

Rate  constant 

(K) 

Solubility 
product  (K£) 

Molar  volume 

(Vk(m3  mol”1)) 

S8(s) 

1.0  s— 1 

19.0  mol  m  3 

1.239  x  10"4 

Li2S8(s) 

1  x  10"11  m6 
mol2  S_1 

18,3400  mol3  m"3 

1.361  x  10"4 

Li2S4(s) 

9.98  x  10"12  m6 
mol2  S_1 

21,480  mol3  m-3 

7.415  x  10  5 

Li2S2(s) 

9.98  x  10"9  m6 
mol2  S-1 

0.6006  mol3  m-3 

4.317  x  10  5 

Li2S(S) 

6.875  x  10  5  m6 
mol2  S-1 

9.95  x  10"4 
mol3  m  3 

2.768  x  10  5 
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